Refining the Spin Hamiltonian in the Spin-i Kagome Lattice Antiferromagnet 

ZnCu3(OH)6Cl2 using Single Crystals 
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We report thermodynamic measurements of the kagome lattice antiferromagnet 

ZnCu3(OH)6Cl2, a promising candidate system with a spin-liquid ground state. Using single crystal 
samples, the magnetic susceptibility both perpendicular and parallel to the kagome plane has been 
measured. A small, temperature-dependent anisotropy has been observed, where Xz/Xv > 1 at high 
temperatures and x^lXv < 1 ^ temperatures. Fits of the high-temperature data to a Curie- 
Weiss model also reveal an anisotropy. By comparing with theoretical calculations, the presence of 
a small easy-axis exchange anisotropy can be deduced as the primary perturbation to the dominant 
Heisenberg nearest neighbor interaction. These results have great bearing on the interpretation of 
theoretical calculations based on the kagome Heisenberg antiferromagnet model to the experiments 
on ZnCu3(OH)6Cl2. 

PACS numbers: 75.30.Gw 75.40.Cx 75.10.Kt 75.50.Ee 
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The quantum spin liquid, a fundamentally new state 
of matter whose ground state does not break conven- 
tional symmetries, has generated much interest in con- 
densed matter physicsii^. It has long been realized that 
the S=^ Heisenberg antiferromagnet on the kagome lat- 
tice (composed of corner sharing triangles) is an ideal 
system in which to look for spin-liquid physics due to 
the high degree of frustration, small spin, and low di- 
mensionality. Herbertsmithite, the x = 1 end member 
of the family Zn paratacamite [Zna;Cu4-a;(OH)6Cl2], is 
arguably one of the best candidate systems to study 
quantum spin liquidsi^ With weak interplane coupling, 
it consists of kagome planes of Cu^"*" ions separated by 
layers of nonmagnetic Zn^+ ions, depicted in Fig. [TJa). 
The current experimental evidence is consistent with the 
presence of a spin-Uquid ground state in this material^^— . 
The Hamiltonian of herbertsmithite consists of a Heisen- 
berg exchange term, with possible perturbations such 
as a Dzyaloshinskii-Moriya (DM) interaction^rJi and ex- 
change anisotropjii^. With a Cu-O-Cu antiferromagnetic 
superexhange interaction of approximately 17 meV, no 
magnetic transition or long range ordering has been ob- 
served down to r = 50 mK— i^ii^. It is important to 
perform measurements on single crystal samples so that 
comparisons can be made to theoretical calculations as- 
suming different perturbations to the Hamiltonia n^"' ^ ^ , 
such as a DM interaction and exchange anisotropy, to 
determine the presence and magnitude of such pertur- 
bations. Resolving this issue is all the more pressing in 
light of recent theoretical work which strongly points to 
a spin-liquid ground state for the S=l/2 kagome lattice 
with isotropic (Heisenberg) exchange- 
Recently, large single crystal samples of the parat- 
acamite family, including herbertsmithite, have been 
synthesizedi^. A powder sample was first synthesized 
inside a sealed quartz tubing and transported under a 
temperature gradient in a three zone furnace for crys- 
tallization. The high quality of the crystals was con- 
firmed by inductively coupled plasma metal analysis, x- 
ray diffraction, neutron diffraction, polarized optics and 
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FIG. 1: (color online) (a) Structure of herbertsmithite with 
Cu^^(big brown spheres) and Zn'^"'" (small red spheres) dis- 
played. The Cu-Cu bonds (thick black solid lines) are all 
equivalent as are the Cu-Zn bonds (thin green dotted lines), 
(b) The oriented single crystal sample (mass — 55.5 mg) of 
herbertsmithite used in the magnetic susceptibility measure- 
ments. 



thermodynamic measurements. Anomalous synchrotron 
x-ray diffraction confirmed the absence of antisite dis- 
order where Zn^"*" ions appear on the Cu sites in the 
kagome layer—. Rather, the main source of disorder is 
the presence of a small fraction of excess Cu^+ ions within 
the Zn inter layers. Raman spectroscopy provides further 
support of a gapless spin-liquid ground stated while /xSR 
points to an easy-axis anisotropy parallel to the c axis 
for magnetizatioE^ii. In this Letter, the magnetic sus- 
ceptibility and specific heat have been investigated with 
fields applied both within and normal to the kagome 
plane. The roles of an easy-axis exchange anisotropy, 
a Dzyaloshinskii-Moriya interaction, and an anisotropic 
g factor for the Cu magnetic moment are discussed. 

Magnetic susceptibility measurements were performed 
on a 55.5 mg single crystal sample of herbertsmithite (2.3 
mm X 2.5 mm x 2.7 mm), shown in Fig. (Hb), using a 
SQUID magnetometer (Quantum Design). The nearly 
cubic shape of the sample minimizes a demagnetization 
correction to the measurements, allowing for a clean mea- 
surement of the intrinsic anisotropy of the material. The 
crystalline axes and the narrow mosaic of the sample were 
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confirmed using an x-ray difi'ractomctcr equipped with an 
area detector. A plastic holder was designed and made 
for securing the crystal for susceptibility measurements 
with magnetic field applied perpendicular (xp) or paral- 
lel ixz) to the crystalline c axis. The background from 
the plastic holder was measured to be negligibly small 
relative to the signal from the sample. In Fig. ^a), the 
quantities XzT and XpT for temperatures between 2 and 
330 K are plotted (where we assume x = M/H in the 
linear regime). The quantity XpowderT, measured on a 
polycrystalline collection of several dozen random orien- 
tated crystals from the same batch, is plotted along with 
XaverageT = ^iXz + '^Xp)T , the Calculated powder av- 
erage. The latter two collapse onto the same curve as 
expected, pointing to the reliability of the single crystal 
measurements. 

In Fig. [2l^b), the anisotropy ratio of the magnetic sus- 
ceptibility calculated as XzlXp is plotted. As tempera- 
ture is increased from 2 to 330 K, the ratio increases from 
0.96 to 1.12 monotonically. The presence of anisotropy 
in the susceptibility agrees qualitatively with susceptibil- 
ity measurements on aligned powdersi^ and recent /iSR 
measurements on single crystalsii. In Fig. [2Jc) , magne- 
tization measurements taken at T = 5 and 300 K are 
plotted as a function of applied field. At T = 5K, the 
anisotropy ratio is close to unity and the two curves over- 
lap for the entire field range. At T = 300 K, there is a 
clear anisotropy with the c axis being the higher suscep- 
tibility direction. The observed magnetic anisotropy is 
independent of the applied field. 

The high quality of the susceptibility data allows for 
further analysis to better understand the intrinsic be- 
havior of the interacting spins on the kagome layers. 
The primary results of this paper are shown in Fig. 3. 
For the susceptibility data at high temperatures, Curie- 
Weiss fits were performed for each data set taken at var- 
ious fields. The Curie- Weiss temperatures and g factors 
determined from the fits (which take into account the 
corrections based on high-temperature series expansion 
calculation o^^i^° ) are plotted in Fig. |3Ja). For both field 
orientations, the Curie- Weiss temperatures and g factors 
increase slightly upon lowering the applied field below 
/io-ff — 0.2T. At//o-ff = IT, the anisotropy ra- 
tio for the g-factor is gz/gp = 1.10. A similar g fac- 
tor anisotropy, though slightly smaller, has been deduced 
from ESR workSi on powders. 

It is important to separate out the anisotropy of the 
Cu moments intrinsic to the kagome planes from that 
related to the impurity spins. The experimentally mea- 
sured magnetic susceptibility originates from both the 
kagome plane and the weakly interacting Cu^^ impuri- 
ties on the interlayer sites. Assuming that the intrinsic 
kagome susceptibility becomes much smaller than the im- 
purity contribution as T — >■ K, consistent with recent 
NMR measurement on single crystal samplea^S, we model 
the impurity susceptibility with a Curie- Weiss law where 
Qcw — 1-3 K^. The best fit gives an estimated 10% 
Cu^+ ions which occupy the interlayer sites for this non- 
deuterated sample. Then, by assuming a temperature 
independent anisotropy for the impurities, the impurity 
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FIG. 2: (color online) (a) Magnetic susceptibility, plotted as 
MT / H = xT, measured under an applied field of /^o-ff = 1 T 
which was oriented both perpendicular to and along the c- 
axis. The data from a powder sample are also plotted and 
compared to the "average" of the single crystal results, as de- 
scribed in the text, (b) The anisotropy x^/Xv of the measured 
susceptibihty plotted as a function of temperature, (c) Mag- 
netization versus field measurements at T = 5 K and 300 K. 
The vertical scale for each temperature is indicated by the 
arrow. 

contribution to the susceptibility can be subtracted re- 
vealing the anisotropy of the intrinsic kagome spins. The 
only remaining free parameter is the anisotropy ratio for 
the impurities (Xz/Xp)imp, and in our analysis, we use 
the value {xzlXp)imp = 1- 

The deduced anisotropy of the susceptibility for the 
intrinsic kagome spins is plotted in Fig. [HJb). The main 
observation, which is relatively insensitive to the model 
parameters, is that XzlXp for the intrinsic susceptibility 
is a monotonically increasing function of temperature for 
T > 150 K. This provides useful information on the im- 
portance of additional terms in the spin Hamiltonian, as 
we discuss further below. Moreover, since we have de- 
duced the anisotropy of the g factor resulting from the 
Curie- Weiss analysis, we can correct for this in deter- 
mining the XzlXp ratio for the intrinsic kagome spins. 
The g- factor corrected data are also plotted in Fig. 3(b). 
At low temperatures (below ~ 5 K) where the impurity 
contribution dominates the susceptibility, the measured 
ratio for Xz/Xp is actually less than 1. If we assume a 
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FIG. 3: (color online) (a) Curie- Weiss temperatures and g 
factors calculated from fits of the magnetic susceptibilities 
between T = 150 and 330 K, as described in the text. The 
proper vertical scale for each data set is indicated by the ar- 
row, (b) The susceptibility anisotropy ratio of the intrin- 
sic kagome spins after subtracting the impurity contribution, 
with and without a correction for the g-factor anisotropy, as 
described in the text. In the model for subtracting the impu- 
rity contribution, temperature independent anisotropy ratios 
{Xz/Xp)imp were assumed. The three curves represent exact 
diagonalization calculationaii for the anisotropy ratio consid- 
ering the effects of an easy-axis exchange anisotropy and a 
DM interaction separately. 

value of {Xz/Xp)iTnp = 0.95, the deduced anisotropy ratio 
for the intrinsic kagome spins exhibits a slight upturn as 
the temperature is cooled below T « 100 K. However, 
the main conclusions of our analysis based on the data 
for temperature above T = 150 K are not quantitatively 
changed. 

Our experimental results shed light on the roles played 
by various perturbations to the spin Hamiltonian beyond 
the nearest neighbor Heisenberg model. The observed 
anisotropy of the intrinsic susceptibility can be compared 
with theoretical calculations using 15-site exact diagonal- 
ization (ED) by Rigol and co-workersii. If an easy-axis 
exchange anisotropy HEA=^^i,jiSf SJ + SfS^) with 
A < is considered, the shape of the anisotropy ver- 
sus temperature curve matches our measurements over a 
wide temperature range, as shown in Fig. 3(b). In fact, 
comparing our ^-factor corrected data with the calcula- 
tion with A = —0.1 J gives a good match for the slope for 
T > 150 K as well as the magnitude for Xz/Xp- The pres- 
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FIG. 4: (color online) (a) Low temperature specific heat data 
on a single crystal sample of herbertsmithite measured un- 
der various applied fields with two orientations, (b) The 
anisotropy ratio of the specific heat measured in the two field 
orientations. 

ence of an anisotropic exchange is consistent with recent 
/iSR measurements on single crystal herbcrtsmithitei^ 
and work on partially aligned powdersi^. The difference 
in the Curie- Weiss temperatures for the two field orien- 
tations is also consistent with the deduced magnitude of 
the easy-axis exchange anisotropy. That is, in Fig. 3 (a), 
Qcw for the field along the c axis is larger than that 
for the field within the kagome plane by about 10%, as 
one would expect for an easy-axis exchange anisotropy of 
A w -0.1 J. 

The DM interaction, HDAi^^t,jD:,{Si x Sj)^ + Dp ■ 
(Si X Sj), has a much smaller effect on the anisotropy 
ratioii. For a wide range of and Dp values (the out- 
of-plane and the in-plane components of the DM vector, 
respectively), the primary effect is to slightly increase the 
anisotropy ratio from unity, where Xz/Xp monotonically 
decreases as temperature increases. The results of two 
model calculations which only include a DM term are 
plotted (one with D^ ~ 0.2 J and one with Dp = 0.2 J). 
Our data appear to rule out such scenarios where only 
a DM term is present, as a small easy- axis exchange 
anisotropy is needed to give the observed Xz/Xp < 1 as 
well as the temperature dependence. 

The specific heat was measured on a 4.10 mg sin- 
gle crystal sample of herbertsmithite using a Quantum 
Design physical property measurement system (PPMS). 
The sample was prepared so that its orientation could 
be changed in situ without remeasuring the background. 
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The specific heat was measured in two field orientations: 
with the field oriented in the kagome plane Cp and per- 
pendicular to the plane Cz ■ The data for a wide range of 
applied fields up to ij^qH = 14 T are plotted in Fig. IH^a). 
The ratio Cz/Cp is plotted in Fig. IH^b) which reveals a 
small magnetocaloric anisotropy. As a check of system- 
atic errors in the measurement, the data measured in 
zero field, taken with both crystal orientations, collapse 
onto the same line. At temperatures higher than 15 K 
(not shown) , the specific heat under all applied magnetic 
fields is identical for both orientations, within error. 

At the lowest temperatures (below ~ 10 K), it is likely 
that the magnetocaloric anisotropy derives from the im- 
purities. The observed Cz/Cp > 1 for ksT ^ fioH indi- 
cates that the impurities are easier to polarize with an in- 
planc field. This idea is further supported by the observa- 
tion in Fig. 2(b) that the ratio Xz/Xp begins to decrease 
very rapidly upon cooling below T = 5 K. This indicates 
the impurity moments develop a (7-factor anisotropy with 
{Xz/Xp)imp < 1 at low temperatures. As another possi- 
bility, it has been shown that the presence of a DM inter- 
action can mix the triplet and singlet stateaiS so that the 
total spin is not a good quantum number. The observed 
field independence of the anisotropy for the susceptibil- 
ity coupled with the field dependence of anisotropy of 
the specific heat point to the possibility that the singlet 
states may be coupled to the applied field. Further the- 
oretical calculations would be useful to determine how 



the thermodynamic quantities should behave under the 
application of in-plane and out-of-plane fields. 

In summary, we have measured the anisotropy of 
the magnetic susceptibility in a single crystal sample 
of herbertsmithite. The temperature dependence of the 
anisotropy allows one to deduce the important additional 
terms in the spin Hamiltonian beyond nearest neighbor 
Heisenberg exchange. A comparison with previous ex- 
act diagonalization calculations indicates the presence of 
an easy-axis exchange anisotropy with A « —0.1 J. This 
small value for the anisotropy indicates that the Heisen- 
berg model is a reasonable approximation to define the 
physics of herbertsmithite. However, calculations start- 
ing from the Ising limit and approaching the Heisenberg 
limit for the S=l/2 kagome antiferromagnet may provide 
useful insight into the behavior of herbertsmithite. A 
field- and temperature-dependent anisotropy in the spe- 
cific heat measured in different field orientations is also 
observed. Further theoretical calculations which include 
an anisotropic exchange interaction as well as a DM in- 
teraction would be most useful for a detailed comparison 
with the data. 
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